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Abstract. As passive optical networks (PON) are increasingly deployed to 
provide high speed Internet access, it is important to understand their funda- 
mental traffic capacity limits. The paper discusses performance models appli- 
cable to wavelength division multiplexing (WDM) EPONs and GPONs under 
the assumption that users access the fibre via optical network units equipped 
with tunable transmitters. The considered stochastic models are based on mul- 
tiserver polling systems for which explicit analytical results are not known. A 
large system asymptotic, mean-field approximation, is used to derive closed 
form solutions of these complex systems. Convergence of the mean field dy- 
namics is proved in the case of a simple network configuration. Simulation 
results show that, for a realistic sized PON, the mean field approximation is 
accurate. 



1. Introduction 

Figure [I] illustrates the main components of a passive optical network (PON). 
Optical network units (ONUs) are situated in or close to user premises. An ONU 
might be dedicated to an individual user or, when fiber is terminated at the curb or 
at the building, it could concentrate the traffic of several users. The ONUs are con- 
trolled from an optical line termination (OLT) equipment that realizes the interface 
with the wide area network. Passive splitters are used to broadcast downstream 
optical signals to the ONUs and to merge upstream signals destined to the OLT. 

To increase the capacity of a PON, it is proposed to employ wavelength division 
multiplexing (WDM) to multiply the capacity of the physical fibre. Typically the 
ONUs will not be equipped to use all wavelengths simultaneously. In this paper we 
assume an ONU has one or more transmitters that can be tuned to any wavelength. 
The spectrum is thus shared efficiently but the capacity of an ONU is limited by 
its quota of tunable transmitter^] 

The major PON design concern is the dynamic bandwidth allocation (DBA) 
algorithm implemented by the OLT to arbitrate use of upstream capacity and to 
avoid collisions. There are two PON standards, EPON from IEEE and GPON from 
ITU, that differ in the way they manage resource sharing. EPON relies on asyn- 
chronous transmission of Ethernet packets while GPON implements a synchronous 
125/xs frame enabling both circuit and packet switching. We develop generic models 
for packet based services that are applicable for both standards. 

The focus of the paper is on traffic capacity, the limit on demand beyond which 
some ONU queues would always be saturated. This is arguably the most critical 
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Figure 1. PON components: ONUs, passive splitters and OLT 



performance for the PON whose very high bit rate ensures delays are tiny at normal 
loads and only become significant as demand attains the capacity limit (see [T1ll3j). 

We assume the OLT manages a single queue per user, receiving queue size reports 
from the ONUs and attributing corresponding grants for upstream transmission. A 
grant is the time interval during which a user may transmit. We maintain that 
such simple management at PON level (i.e., no accounting for intra-ONU class 
of service differentiation at the OLT) is highly desirable to limit complexity and 
generally meets user requirements, the OLT ensures users receive a fair share of 
bandwidth, the ONUs realize what service differentiation is needed within that 
share by choosing which packets to send to fulfill the allocated grants. 

The notion of fair allocation naturally leads to a form of polling. For EPON we 
define a DBA that realizes a classical polling system with a particular interpretation 
for the switchover time between queues. The same model can be applied to one 
possible DBA implementation for GPON. A second considered GPON DBA, on the 
other hand, is more naturally modelled as a head of line processor sharing system. 
See below for more detail. 

The underlying polling systems are complicated by the fact that ONU capacity 
is limited by the number of tunable transmitters they have, i.e., the number of 
servers that can be present simultaneously at an ONU. In particular, stability con- 
ditions of multiserver polling systems with limited-gated policies appear intractable 
in general. 

Related work. There is an abundant literature on PONs and numerous DBA al- 
gorithms have been proposed. Two recent surveys, by McGarry et al. [16] and 
by Zheng and Mouftah [23], present a comprehensive review for EPON. Skubic et 
al. [IH] compare DBAs used in EPON with alternatives proposed for GPON. The 
DBA algorithm we propose is similar in principle to IPACT, proposed by Kramer et 
al. [TU] for a regular EPON and extended by Kwong et al. for WDM [IT]. More 
complex class of service aware DBAs for a WDM EPON are proposed in [15l [T4J 
[131E1II7]. A paper by Gliwan et al. discusses DBA for a coarse WDM GPON [8]. 
Song et al. recognize that classical EPON DBAs can lead to lost capacity in long 
reach EPONs and suggest the ONUs should initiate multiple interleaved cycles or 
threads [2U] . We advocate an alternative DBA that avoids loss of capacity whatever 
the propagation times [TJ. 

Most performance evaluations of DBA algorithms rely on simulation. A few 
authors have developed analytical models, generally using techniques developed for 
polling systems. Park et al. derived closed form formulas for the mean packet delay 
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in a symmetric PON with identical ONUs and negligible propagation times under 
the assumption of Poisson traffic and gated service discipline [T8]. The two-stage 
polling system identified in that paper has been further analyzed by van der Mei 
and Resing [22] when ONUs have heterogeneous load. To the best of our knowledge, 
polling models have not been developed explicitly for WDM PONs. The capacity 
of classical multiserver polling models was analysed by Fricker and Jaibi [7 . Borst 
and van der Mei introduced the notion of server limit in [3] but did not evaluate 
capacity limits. Down [5] presents several capacity results for multiserver polling 
models with server limits but with unlimited service policies. Lastly, the study 
of head of line processor sharing systems by Brandt and Brandt is relevant but 
unfortunately provides no exact capacity results [4]. 

Mean Field Limits of Polling Systems: To overcome the complexity of the polling 
systems associated with these networks, the main contribution of the paper consists 
in the analysis of a large system limit where both the number of ONUs and the 
number of wavelengths grow indefinitely together. This scaling, also called mean 
field limit, provides explicit formulas and simple numerical algorithms to derive 
useful approximations for realistically sized PONs. It is based on the fact that, 
when the approximation holds, the ONUs become statically independent in the 
limit, their interaction being represented through a deterministic equation. See |21j 
for a general presentation of these methods. Mean field limits have been used for 
some time now in queueing theory notably to study circuit switched networks, 
see [S] for example. 

In Section [2] we outline our proposed DBAs for EPON and GPON. We then 
recall relevant multiserver polling system results and discuss the difficulty of an 
exact analysis of the WDM PON. The large system asymptotic is developed in 
Section [4] where practically useful traffic capacity results are derived using the mean 
field assumption. This property is proved for a particular network configuration in 
Section [5j The validity of this assumption is tested, in penultimate Section [6] and 
used to evaluate the impact of switch overhead on PON traffic capacity. Section [7] 
concludes the paper. 

2. Dynamic bandwidth allocation 

The DBAs we consider for both EPON and GPON relate to a single, classless 
queue per user. The objective is to ensure low latency by imposing grant limits per 
visit and to ensure users receive adequate throughput by controlling fairness. 

2.1. EPON DBA. The EPON DBA is governed by the multi-point control pro- 
tocol (MPCP) based on the exchange of REPORT and GATE messages between 
OLT and ONUs, as specified in IEEE 802. 3ah. An extension to WDM is proposed 
by McGarry et al. [T5] . ONUs report current contents of their user queues via RE- 
PORT messages while GATE messages are sent downstream by the OLT to allocate 
grants for each queue. In [T], we proposed a novel DBA where the OLT precisely 
times the sending of GATE messages and the upstream return of data packets and 
REPORTs to maximize traffic capacity. The timing is designed to compensate for 
differences in propagation times and to avoid wasting capacity when propagation 
times are large. 

The EPON behaves then like a polling system. The ONUs are "visited" by 
each wavelength following a certain schedule. The visit enables the transmission 
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of upstream data and a REPORT. The grants account for maximum per queue 
quotas. Each visit incurs a switch overhead time equal to the time needed to send 
the REPORT plus a physical layer guard time. The overhead is around 2/xs for 
a lGb/s EPON [IS]. The OLT is aware of all allocations and readily implements 
constraints on the number of simultaneous transmissions an ONU can sustain. 

2.2. GPON DBA. The GPON standards define a particular framework for DBA 
that allows resource sharing between several communication services (see the G.984 
series of ITU recommendations). We consider only the packet switched service. 
Ethernet packets are transmitted within GPON transmission convergence (GTC) 
layer frames of 125/xs fixed duration. Downstream GTC frames are also used to 
convey grants from OLT to ONU users and reports on queue contents are included 
in upstream GTC frames from ONU to OLT. Grants in a downstream GTC frame 
attribute transmission time slots in a precise upstream GCT frame occurring after 
a fixed time offset. The offset is calculated to allow all ONUs to fulfill their grants 
accounting for processing time and propagation delays. 

It is possible to design a DBA in this framework that emulates a polling system. 
Wavelength capacity is allocated sequentially to ONUs for data transmission, or just 
for reporting if the queues are currently empty. Overhead is incurred on switching 
from one ONU to the next, as for the EPON. However, this overhead is rather 
small according to the GPON specification, accounting for only some lOOrys [T§] . 
An additional fixed overhead is incurred at the start of each frame. 

The GPON DBA might alternatively consist in providing every ONU queue a 
share of every upstream GTC frame on one or several wavelengths. The OLT 
would calculate an allocation for each queue and a schedule of transmissions based 
on previously received reports and accounting for ONU transmission capacity. In 
this way, there is a fixed overhead per frame, accounting for reports and physical 
constraints, to be subtracted from the overall capacity. The remaining capacity 
would be shared between users according to some fairness policy. 

3. MULTISERVER POLLING SYSTEMS 

We introduce notation and discuss known polling system results that can be used 
to analyse WDM PONs. 

3.1. Notation and assumptions. In the following wavelengths are identified with 
servers circulating among ONUs. The polling system has N ONUs and L wave- 
lengths. ONU i is equipped with U tunable transmitters and manages ni indepen- 
dent user queues, i.e., servers can be present at ONU i simultaneously. When 
the DBA emulates a polling system, the j th queue of ONU i, queue has a 

maximum grant of dij seconds per visit by any wavelength. This mechanism is 
known as the limited-gated service policy. Finally, Aj is the switch overhead time 
associated with ONU i. 

We make no assumptions about the traffic in the user queues except that the 
packet arrival process is stationary and has a well defined rate. The traffic intensity 
of queue is ptj expressed in units of the bit rate of one wavelength. We define 

Pi =J2jPij and P = Y.i Pi- 
Pi. wavelength is assumed to visit ONUs in a fixed cyclic order, assumed differ- 
ent for each wavelength. For periodic polling, visits to ONUs occur according to a 
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deterministic schedule. A wavelength that should visit an ONU with no free trans- 
mitter is assumed to skip directly to the next ONU in the schedule. Alternatively, 
for random polling, the next ONU to visit is chosen at random from the ONUs that 
currently have a free transmitter. 

3.2. Limited grants, no server limits. In this section we assume U = L for all i 
- an ONU can be served by an arbitrary number < L of servers with a limited-gated 
policy. Polling is either periodic or random. 

Proposition 1. All queues in all ONUs are stable if and only if, for all i,j, 

N 

(1) p+^^A,<i. 

a v k=i 

When not all queues satisfy |Ip, there is a subset S whose members are saturated 
while, for 4- we have 

(2) 9 - E ^ + S(f> + £ d A< L - 

(Lk)eS lj \k = l (Lk)eS J 

Proof. Condition derives from Theorem 3 in [7j or Theorem 2.4 in [6j with 
appropriate changes in notation. Relaxed conditions in ([2| are deduced from the 
discussion in Section 4 of [7]. In fact, the results of [7] and [6] must be modified 
slightly since grant limits here apply to transmission time and not packets but this 
is straightforward. □ 

We later require the following additional result. 

Proposition 2. When the system is stable, the mean time C between successive 
visits of wavelengths to any given ONU is C — (Ai + A2 + • • • + Ajy)/(L — p). 

Proof. The result may be deduced for periodic polling from Section 3 in 3 . The 
following alternative proof is valid also for random polling. Since the system is 
stable, the mean number of wavelengths transmitting data is p. Thus L — p is the 
mean number of wavelengths currently busy due to switch overhead. The fraction 
of time a server at ONU i experiences a switch overhead is Aj/C. The relation for 
C is obtained by summing up these quantities. □ 

3.3. Server limits, no grant limits. The capacity of a multiserver polling system 
with server limits but unlimited gated service with periodic or random polling is 
given by the following. 

Proposition 3. When the WDM EPON with tunable transmitters has unlimited 
gated service, it is stable if pi < ti for 1 < i < N and p\ + pi + • • • + pN < L. 

Proof. The conditions are clearly necessary since otherwise at least one user queue 
would be unstable. To prove sufficiency, we use fluid limit arguments as in 
Starting from any initial conditions, it is necessary to prove that the overall fluid 
content denoted W(t) decreases with time. Let J"(£) designate the number of ONU 
queues that have a fluid backlog at t. If X^eJ^t) ^ — L, all the service capacity 

is devoted to fluid queues and W(t) = J2i Pi ~ L < by the second condition. If 
Siej^(t) *t < L, all queues in F(t) will be served at their maximum rate. We have, 

therefore, W(t) = E <6 ^(t)(« ~ **) < °- D 
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Figure 2. Sharing L wavelengths between Z(t) active queues - no 
switch times 

Combining server limits and grant limits leads to a system whose traffic capacity 
turns out to be intractable, even for a simple toy example with 3 ONUs with a 
single queue, 2 wavelengths and no switch overhead pQ. This motivates the large 
system approximations discussed in the next section. 

4. Large system asymptotic 

In this section wc derive the traffic capacity of a WDM PON with tunable 
transmitters under the assumption that a mean held limit is valid in large systems. 

4.1. Negligible switch overhead. As a preliminary, consider the polling system 
representing a WDM PON under the assumption that overhead (for switching be- 
tween ONUs) is negligibly small. Assume the ONUs are each equipped with one 
tunable transmitter. Let the number of active ONUs at time t be Z(t), i.e., Z(t) is 
the number of ONUs with at least one backlogged queue. When Z(t) < L, every 
ONU is served by one wavelength. If Z(t) > L, on the other hand, the PON shares 
the capacity of the L wavelengths according to the DBA algorithm. 

Figure[2]illustrates the importance of the value of L. In the top figure, L = 3 and, 
depending on the ONU loads, it is a frequent occurrence that Z(t) > L. ONUs then 
share capacity dynamically and it is not easy to predict PON performance. In the 
lower figure, L = 15 and the fluctuations of Z(t) have relatively lower amplitude. 
In the depicted realization, every PON is always able to use its tunable transmitter 
whenever it is active. 

As L and N increase indefinitely, it is intuitively clear that the fluctuations of 
Z(t) become small compared to L. Each ONU then behaves like an independent 
polling system with the server capacity of one wavelength. Every ONU queue is 
stable if p t < 1 for 1 < i < N and p < L. 

This intuition is confirmed by the propositions of Section [5] in the case where 
ONUs have a single queue, service is limited to one packet per cycle, packet arrivals 
are Poisson and packet sizes exponential. A more general proof is, for the moment, 
out of reach. We nevertheless proceed in this and the next sections with the rea- 
sonable assumption that a "mean field" asymptotic is valid as L and N both tend 
to infinity. 
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Figure 4. ONU local polling system 



It is clear that the above capacity limits can be generalized to ONUs equipped 
with a variable number of tunable transmitters. In a large system, each ONU always 
has access to its quota of wavelengths whenever it is active. The traffic capacity 
of ONU i is thus that of a ^-server polling system: queues are stable if pi < t{ for 
1 < i < N and p < L. 

4.2. Limited service polling system with switch overhead. We now suppose 
switch overheads are non-negligible: each visit to ONU i consumes Aj seconds 
of wavelength capacity. Again, first assume each ONU is equipped with just one 
tunable transmitter. Figure [3] shows how wavelengths are fully used in a WDM 
PON. Each wavelength is either transmitting data or is unavailable because of the 
switch overhead. When p < L, the visit frequency to inactive ONUs is such that 
residual capacity L — p is entirely consumed by overhead. 

Figure [4] illustrates the upstream activity of an ONU. This ONU has 3 queues 
that are served cyclically respecting maximum grant limits on each visit. In the 
figure, queue 3 happens to be empty in cycle n. Each set of grants terminates with 
the switch overhead. After the overhead, the wavelength visits another free ONU 
(possibly this one again) and the considered ONU sees a vacation time (> 0) before 
it is again visited by a wavelength. 

When the system is stable, at equilibrium the vacation time of ONU i, i.e., the 
time when no server is present, is denoted by O^. Note that this random variable 
depends a priori on the activity of all ONUs. Under the mean field assumption 
below, this dependence vanishes in the large system limit where N and L tend to 
infinity. Let 9 = E[0,-] and 5 be the switchover time at equilibrium over all ONUs. 
First, we derive a useful general relation. 



Proposition 4. At equilibrium, the mean vacation interval i 
tunable transmitter per ONU satisfies 

n (N-L)S 



of a system with one 



(3) 



P 



Proof. At equilibrium, the number of ONUs that are not currently being served 
is constant and equal to N — L. Consequently (N — L)/8 is the rate at which 
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one of these ONUs receives the visit of a server. Similarly, the expected number of 
wavelengths transmitting data must be equal to the offered traffic p. The remainder, 
L — p, is the mean number of servers in overhead period and consequently (L — p)/5 
is the rate at which servers leave a switch overhead period. The identity of the 
proposition follows. □ 

To exploit Equation ^ we must make the following additional assumption. 
Mean field assumption We assume successive intervals 0™ are independent and 
identically distributed for all i and all n with mean 9. This means the local ONU 
polling system with vacation depicted in Figure|4]can be analyzed as an independent 
system, the impact of the other ONUs being manifested through the value of the 
mean field parameter 9. 

We maintain this assumption is reasonable when N and L are large and when, for 
each wavelength, the choice of the next ONU to visit is made uniformly at random 
among those that are eligible. The latter assumption applies for both random 
polling and periodic polling in the large system limit. This, roughly speaking, is 
because the past activity of any ONU has negligible impact on future wavelength 
visits while available destinations for a wavelength's next visit hardly depend on its 
own history. The system is effectively memoryless and, indeed, we would further 
claim the distribution of the 0" is exponential in the large system limit. This 
property is not used for the present capacity results. 

To evaluate 9, note that, under the mean field assumption, the system repre- 
sented in Figure [4] behaves like a classical single server periodic polling system with 
limited gated service and expected overhead per cycle equal to A, +9. By Proposi- 
tion [l] its mean cycle time is (Aj + 9) / (1 — pi) . Straightforward calculations show 
that the overall mean overhead S defined above can then be expressed as 

m c = E(l-Pi)A«/(A i + 0) 

1 ) E(i-ft)/(A, + ^ ■ 

Eliminating S in ^ and Q yields an equation that can be solved for 9. When 
the Aj are all equal (to 5), 9 is given directly by ([3|. Otherwise, we must resort to 
numerical evaluation. 

Stability conditions result from those of the individual ONU polling systems 
where the overhead and vacation are assimilated to the per cycle switch time. 

Proposition 5. Under the mean field assumption, the WDM PON with one tunable 
transmitter per ONU, represented as a multiserver, limited-gated polling system 
where no more than one server can simultaneously attend any queue, is stable if 

(5) Pi + ^-(Ai + 6) < 1, l<i<N, l<j< ni . 

where 9 is derived from |5|) and Q). 

The proof follows directly from Proposition [T] on making the appropriate sub- 
stitutions. When the overhead is the same for each ONU we readily deduce the 
following. 

Corollary 1. When Aj = 8 for all i, every queue is stable if 

(6) p + max(g&^)<2, 



'j 
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4.3. Local stability. Limited gated service is designed to limit the impact on 
ONU performance when some users saturate their share of PON capacity. Let the 
set of saturated queues be S. To derive stability conditions for the other queues we 
must simply replace Aj by Aj + j)eS m Q to compute 9 and substitute 
for Ai and 6 in §5§. 

For a system with given traffic distribution, the following algorithm determines 
the set S of saturated queues. 

(1) S = 4>, the empty set. 

(2) find = argmax (lj) ^ 5 {giAJ. 

(3) if ^ is satisfied for then stop. 

(4) else, 

- Ai <- A % +dij, 

- 9 <- P~ Pij, 

- S<-SU(i,j). 

(5) repeat from step[2j 

4.4. Multiple transmitters. To simplify the presentation, we have so far assumed 
each ONU has exactly one tunable transmitter. However, the above results apply 
immediately with slight modification to a variable number of transmitters if we 
make the following assumption. The wavelengths on terminating one service visit 
an eligible ONU with probability proportional to its number of free transmitters. 
This would result in the large system limit if each wavelength had a fixed cycle 
of visits chosen as a random permutation of transmitters. In condition ([5]), it is 
necessary to replace 1 on the right hand side by i,. 

4.5. Frame based GPON DBA. In GPON, a possible DBA is to divide the 
capacity available in a frame (or several successive frames) between ONU queues, 
applying a max-min fair algorithm. The overhead per frame is fixed and equal to 
J^Aj. The residual capacity would be shared by the OLT between active ONU 
queues. Queue j of ONU i would receive an allocation when active proportional to 
its weight w^. 

Assume the frame duration / and the Aj tend to zero while the ratio Aj// 
remains equal to Si. The underlying service system is now weighted head of line 
processor sharing with a maximum service rate for ONU i equivalent to U(l — Si). 

The capacity of this system appears intractable in general [4] . In the large system 
limit where the number of queues and the system capacity increase indefinitely, 
however, we intuitively expect the service rate constraint to apply to all ONUs, i.e., 
with high probability, every ONU is served when busy at its maximum rate. The 
system in fact behaves as discussed in Section |4.1| and the stability conditions are 

(7) p < L(l - ^2 Si), and pi <ti(l-8i). 

5. Proof of Mean Field Convergence 

The above results rely on the assumption that, as N and L get large, the indi- 
vidual ONU queue sets behave as independent systems. This is a classical mean 
field assumption that would ideally need formal justification. In this section we 
prove the mean field property under some further simplifying assumptions. 
Homogeneous limited server polling system without switch overhead times. We con- 
sider an S Ar-server polling system where each of N queues can be served by only one 
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server at a time. After each service, a server chooses a new, non-empty queue uni- 
formly at random. Moves between queues are instantaneous, there is no switchover 
time. Customer arrivals are Poisson of intensity A in each queue and service times 
are exponential with rate p. The queue load is p = X/fi. 

For t > 0, let Qi(t) be the number of customers in the ith queue at time t and 
h{t) £ {0, 1} an indicator of the presence of a server at queue i. Clearly, if Qi(t) = 
then 2j(t) = 0. The process 

(X(t)) d ^ ((Q n (t),I n (t),l<n<N)) 

is then a Markov process taking values in S = {x = (x n , i n ) e (Nx {0, 1})^ : x n > 
0ifi n = l}. 

We consider a "large" system where 

lim = s, with s G (0, 1). 

JV-H-oo N 

Note that the "natural" necessary stability conditions are Np < Sn or, in the limit, 
p < s. 

First assume the initial occupancies Qi(t) of each queue are independent and 
geometrically distributed with parameter p. As N grows large, we have 

JV 

E 1 {Q»(o)>o> ~ WP(Qi(0) > 0) = Np < Ns ~ S N . 

n=l 

Under this regime there is an infinite number of idle servers so that the vector 
(7„(0)) is uniquely determined. The following proposition shows that this holds 
with high probability on any finite time interval. 

Proposition 6. If p < s and the initial state (X(0)) = ((Q„(0), / n (0)) is such that 
the variables Q„(0) are independent with a geometric distribution with a parameter 
p then, for each T > 0, there exists e > such that 

/ N \ 

,E 



lim P sup Vl {Q „ (t)>0} <SW-£7V =1. 



Proof. A simple coupling argument is used. Consider N independent M/M/l 
queues with lengths L n (t) with arrival rate A and service rate p at equilibrium. 
For each t > 0, L n (t), 1 < i < N are then N independent geometric random 
variables with parameter p. Let 

1 N 

A N(t) = Z2 1 {L n (t)>0}- 
n=l 

By the law of large numbers, almost surely, 

lim Apf(t) = p. 

In order to prove a stronger statement, i.e., that this convergence holds not only 
for a fixed time but on a whole time interval, one has to work in the Skorokhod 
space of probability distributions on real valued functions on [0, T] with limits on 
the left and continuous on the right. See Billingsley [2], for example. 

In the following, (Af x ^) denotes an i.i.d. sequence of Poisson processes with 
parameter x and Af x ,i([a,b]) is the number of points in the interval [a,b]. We first 
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prove that the sequence of processes (A]y(t)) is tight. Since a large variation of 
(A N (t)) can only be due to either arrivals or departures, we have, for 5 and r\ > 0, 

P[ sup \A N (t) - A N (s)\ > n ) 

\l<s,t<T:\t-s\<6 I 



( 1 N 

<P sup t^-AMM) > 



\l<s,t<T:|t-s|<« ^ n=1 -y 

+ P( Sup ^ 5Z^,n([«. *])>?)■ 

\l<s,t<T:\t-s\<S iV ^ 1 j 

Using the fact that A/a,i H + N\,n has the same distribution as AFnx.i, by the 

law of large numbers for Poisson processes, the first term of the right hand side of 
the above expression becomes arbitrarily small as N grows if S is less than tj/2X. 
Similarly, the second term has a similar property. We deduce from this estimate 
that the sequence of processes (Ajv(t),0 < t < T) is tight for the convergence of 
the uniform norm and, therefore, that any of its limiting points is a continuous 
process. Since the marginals converge to the constant p, the only limiting process 
is the constant process equal to p. 

The sequence (Ajv(£),0 < t < T) thus converges as a process to (p). In par- 
ticular, the supremum of (Ajv(t)) on the time interval [0, T] converges to p. Since 
Sn ~ sN > pN, for e < s — p, we have 

lim P ( sup A N (t) < ^ - £ ) = 1. 
JV^+oo Vo<t<T JV / 

On the event {sup 0<t<T Ajv(t) < Sn/N — e}, there is always an idle server. In this 
situation, polling system (X(t)) is equivalent to N independent M/M/l queues so 
that (L n (t),0 <t <T) has the same distribution as {Q n (i),0 <t<T). □ 

The following corollary states the required mean field property given the assumed 
initial occupancy distribution. 

Corollary 2. Under the assumptions of the above proposition, on any finite time 
interval, as N gets large, the Markov process (X(t)) has the same distribution as 
N independent M/M/l queues. In this limit, the natural condition p < s is the 
stability condition of the system. 

We conclude this section by a result which states that, starting from any initial 
state, the process very rapidly enters a set of states where many servers are idle. 

Proposition 7. If p < s and x = (q n ,i n ) € S is such that \{n : q n > 0}| > Sn 

and qi + ■ ■ ■ + qN < CN for some constant C, the infimum 

T N = inf{s, \{n,Q n (t) > 0}| < S N }, 

satisfies 



Proof. Define, for t > 0, 



N xy ~ hSn/N-X' 
Z(t) = Q 1 (t) + --- + Q N (t). 
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Since all servers are busy at time by assumption, as long as t < T/v, Z(t) — Z(0) can 
be represented as the difference of two Poisson processes with respective parameters 
XN and /zSjv, so that (Z(t A T/v) + (/aSn — NX)(t A T/v)) is a super-martingale. In 
particular 

E(Z(t A T N ) + (pS N ~ NX)(t A T N )) < Z(0) < CN. 

Hence, (jiSn — NX)E(t A Tjv)) < CN and the desired inequality follows on letting 
t go to infinity. □ 



6. Evaluation 

We report results of a preliminary evaluation of the accuracy of the mean field 
approximation and quantify the impact of the switch overhead. 

6.1. Convergence to mean field. In the simulations, we consider a realistic PON 
configuration where each wavelength provides an upstream transmission bit rate of 
1 Gbps. We assume that ONUs have a single queue with infinite capacity and one 
tunable transmitter. Users send fixed size 1000 byte packets according to a Poisson 
process. The grant size is fixed to 8 /is (transmission time of one packet) and switch 
overhead is set to 1.2 fis (15% of the grant size). 

Figure [5] shows the proportion of wavelengths in switch overhead measured in 1ms 
intervals for N = 10, 50, 100 and 500, respectively. Overall load is 0.2iV and L/N = 
0.5. The figure shows that, as the number of ONUs increases, the fluctuations of 
this proportion decrease. The wavelength visit rate to a given ONU, determined 
by this proportion, thus converges to a constant. This is why it is plausible that 
individual ONUs tend to become statistically independent. 

6.2. Capacity region. Figure [6] depicts the capacity region of a PON with two 
classes of ONUs. Each ONU of the same class has the same load and there is an 
equal number of ONUs per class. L/N = 1/2. The region to the left and below 
the plots corresponds to loads for which the ONU queues are stable. Units are 
normalized loads with respect to the total number of ONUs. 

Figure [6] compares the capacity region determined by ^ with simulation results 
for TV = 6 and N = 20. The accuracy of the results is good for N — 6 and excellent 
for N = 20. Discrepancies for N = 6 with unbalanced traffic are not surprising 
as the number of ONUs is then effectively only 3. Similar precision is observed in 
Figure [7] where class 2 has three times as many ONUs as class 1. Note that PONs 
are currently designed for around 100 ONUs so that the accuracy of the mean field 
approximation is largely sufficient. 

6.3. Impact of overhead. Figure [8] illustrates the impact of the overhead as a 
percentage of the grant limit given by the mean field approximation. For these 
results, L/N = 0.4 and 40% of ONUs are in class 1. The graphs show that the 
impact of overhead is significant. It is worthwhile increasing the grant size as long 
as the worst case latency remains sufficiently small. Note that when some queues 
are overloaded, their grant is effectively added to the overhead seen by the other 
queues, as discussed in Section [43] 
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FIGURE 5. Proportion of wavelengths in switch overhead in in- 
tervals of 1 ms for different values of ONUs (N = 10, 50, 100 and 
500) 




Normalized load (class 1) 



Figure 6. Stability region for two classes of ONUs with balanced 
number of ONUs 
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Figure 7. Stability region for two classes of ONUs with unbal- 
anced number of ONUs 




Figure 8. Stability region for two classes of ONUs with unbal- 
anced number of ONUs for different ratios switchover/Grant Size 



7. Conclusion 

The design of dynamic bandwidth allocation algorithms for WDM PONs brings 
challenging performance problems, even when we suppose the DBA does not at- 
tempt to realize intra-ONU service differentiation. The main performance criterion 
for these systems is the maximal traffic capacity they can achieve, delays being 
negligibly small until loads are very close to this capacity. 

The natural stochastic models to be investigated in this domain are multiserver 
polling models with both server limits and grant limits. They appear to be not 
tractable in general. For this reason, we have explored a large system approximation 
based on a mean field assumption. Under this assumption, the behaviour of a 
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given ONU can be described as an isolated queueing system independent of the 
other ONUs. The traffic capacity analysis for this limiting representation of the 
PON notably reveals the non-trivial impact of switch overhead. Simulations reveal 
excellent accuracy even for systems much smaller than that of real PONs. 

The mathematical justification of the mean field approximation has been pro- 
vided in a simple case where switch times are null and traffic is Markovian. In 
future work we aim to extend this analysis to the realistic case of general traffic 
and non-zero switch overhead. We also intend to use the results derived here to 
perform a more thorough analysis of DBA performance for EPON and GPON. 

Note finally that the derived large system approximations can be applied in other 
contexts where the natural model is a multiserver polling system with both server 
and grant limits. 
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